Highly Nonlinear Solitary Waves in Periodic Dimer Granular Chains 
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We report the propagation of highly nonlinear solitary waves in heterogeneous, periodic granular 
media using experiments, numerical simulations, and theoretical analysis. We examine periodic ar- 
rangements of particles in experiments in which stiff er/heavier beads (stainless steel) are alternated 
with softer/lighter ones (polytetrafiuoroethylene beads). We find excellent agreement between ex- 
periments and numerics in a model with Hertzian interactions between adjacent beads, which in 
turn agrees very well with a theoretical analysis of the model in the long-wavelength regime that 
we derive for heterogeneous environments and general bead interactions. Our analysis encompasses 
previously-studied examples as special cases and also provides key insights on the influence of the 
dimer lattice on the properties (width and propagation speed) of the obtained highly nonlinear wave 
solutions. 
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Over the past several years, the highly nonlinear dy- 
namic response of granular materials has drawn increased 
attention from the scientific community [l|, [2], y, |J, |5|, |6|, 
0, H, 0, El EI El, El El El EH- The corresponding 
theory developed for uniform lattice systems [l[ supports 
the formation of a novel type of wave in materials, setting 
a paradigm for the design and creation of systems with 
unprecedented properties. A simple setup for the study 
of highly nonlinear dynamics in solids is provided by 
one-dimensional (ID) granular media consisting of chains 
of interacting spherical particles that deform elastically 
when they collide. The broad interest in such systems 
has arisen because they possess qualitatively different 
features from weakly nonlinear systems. For example, 
their solitary-wave solutions have a finite support that is 
independent of their amplitude , providing perhaps the 
most experimentally tractable application of the notion 
of "compactons" [17|. There have also been a number 
of recent studies on the effects of defects (i.e., inhomo- 
geneities, particles with different masses, etc.) in such 
systems, allowing the observation of interesting physical 
responses such as fragmentation, anomalous reflections, 
and energy trapping 
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over, chains of granular media have been shown to be 



ing with the Fermi-Pasta-Ulam problem [18|, [ljj, |2fJ • Its 
applications arise in numerous areas of physics, includ- 
ing coupled waveguide arrays and photorefractive crys- 
tals in nonlinear optics 



21, \22 



Bose-Einstein conden- 
sates in optical lattices in atomic physics [23[ , and DNA 
double-strand dynamics in biophysics [2J ]. A particular 
theme that often arises in this context is that of "het- 
erogeneous" versus "uniform" lattices. Here, we focus on 
the prototypical heterogeneity of "dimers" (i.e., chains 
of beads made of alternating materials). This topic is 
of interest to a diverse array of physical settings rang- 
ing from ferroelectric perovskites [25L 1261 ] and polymers 
[27| to optical waveguides [2§| and cantilever arrays . 
Chains of beads provide a particularly interesting and 
versatile setting with which to study such problems be- 
cause of their strong nonlinearity and the wide range of 
available material properties (and concomitant tunabil- 

ity) HIElEir 



II E|- More- 
over, chains ot granular media nave been shown to be 
highly tunable 0, 0, 0] and have the potential to be used 
in many engineering applications - including shock and 

sound focusing 
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energy absorbing layers 
devices (tunable acoustic lenses and delay lines), sound 
absorption layers, and sound scramblers 0, E3|- 

A class of phenomena for which ID chains of granular 
media provide an ideal setting concerns the interplay be- 
tween nonlinearity and periodicity. The study of nonlin- 
ear oscillator chains has a time-honored history, originat- 



Our aim in the present work is to investigate solitary 
wave propagation in chains of granular dimers [301 ] com- 
bining experiments, numerical simulations, and theoret- 
ical analysis. Some preliminary results on dimer chains 
are discussed in [lj. In this Rapid Communication, we 
consider a broad class of configurations consisting of 
stainless steekPTFE (polytetrafiuoroethylene) [H, E3 
dimer chains with different periodicities (obtained by 
varying the number of consecutive steel particles). We 
report very good agreement between experiments and 
numerics. We also apply a long- wavelength approxima- 
tion to the nonlinear lattice model to obtain a quasi- 
continuum nonlinear partial differential equation (PDE) 
description of the system. Because of the chain lengths 
and the pulse amplitudes considered in our experimen- 
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FIG. 1: (Color online) (a) Experimental setup for dimer chain 
consisting of a periodic array of N consecutive stainless steel 
beads interspersed with 1 PTFE bead, (b) Schematic diagram 
of the composition of the sensors placed in the chain. 
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FIG. 2: (Color online) Force versus time response obtained 
from chains of dimers consisting of (a,b) 1 or (c,d) 2 stain- 
less steel beads alternating with 1 PTFE bead. Panels (a,c) 
show experimental results, and (b,d) show the corresponding 
numerical data. The striker has a velocity of 1.37 m/s on im- 
pact and the y-axis scale is 2 N per division. The numbered 
arrows point to the corresponding particles in the chain. In 
both configurations, the second curve (showing the results for 
particle 24) represents a PTFE bead and the other curves 
represent steel beads. 



tal and numerical analyses, the pulse propagation can be 
safely assumed to be within the highly nonlinear regime 
and the effects of gravity in the theoretical analysis have 
been neglected, in accordance with the discussion of 0. 
We obtain analytical expressions for wave solutions of the 
PDE and find very good agreement between the widths 
and propagation speeds of these solutions with those ob- 
tained from experiments and numerical simulations. 

Experimental Setup. The experimental dimer chains 
were composed of vertically aligned beads in a delrin 
guide that contained slots for sensor connections. Each 
"N : 1 dimer" (see Fig. [lj included a variable number 
N G {1, . . . , 7} of the high- modulus, large mass stainless 
steel beads (non-magnetic, 316 type) alternating with a 
single low-modulus, small mass PTFE bead in a periodic 
sequence. The diameter of all spheres was 4.76 mm, and 
the number of beads in the chain was 38. Three piezo- 
sensors were embedded inside the particles in the system, 
as described in 

US US El ■ The calibrated sensors (RC 
« 10 3 fis) were connected to a four-channel Tektronix 
Oscilloscope (TKTDS 2014), allowing direct visualiza- 
tion of the propagating pulse via force versus time curves 
and time-of-flight calculations of the pulse speed. Waves 
were generated by a striker (a 0.45 g stainless steel bead) 
dropped from various heights. The steel beads had mass 
0.45 g, elastic modulus 193 GPa, and Poisson ratio 0.3 



tween beads [l[: 
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—j ; and the PTFE beads had mass 0.123 g, elastic 
modulus 1.46 GPa, and Poisson ratio 0.46 [H, |H3]. Pre- 
vious studies of dimer chains only considered materials 
with similar elastic moduli 0]. 

Numerical simulations. We model a chain of n spher- 
ical beads as a ID lattice with Hertzian interactions be- 
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where yj is the coordinate of the center of the jth particle, 
j e {1, ■ • • ,n}, Sj ee max{yj-i-yj,0} for j e {2, . . . ,n}, 
Si ee 0, 5 n+ i ee max{y„,0}, g is the gravitational accel- 
eration, Ej is the Young's modulus of the jth bead, Vj 
is its Poisson ratio, rrij is its mass, and Rj is its radius. 
The particle j = 1 represents the steel striker, and the 
(n + l)st particle represents the wall (i.e., an infinite- 
radius particle that cannot be displaced). The initial ve- 
locity of the striker is determined by experiments and all 
other particles start at rest in their equilibrium positions. 

In Fig. [21 we show the very good agreement between 
experimental and numerical results for 1 : 1 [panels (a,b)] 
and 2 : 1 [panels (c,d)] dimers of steeLPTFE particles. 
For the 1 : 1 dimers, the dynamics indicate that the initial 
excited impulse develops into a solitary wave within the 
first 10 particles of the chain. We also obtain robust 
pulses for N : 1 dimers with N > 1 (with N as high as 7 
in experiments and as high as 22 in numerics), though the 
transient dynamics and spatial widths of the developed 
solitary-like waves are different. 

Theoretical Analysis. We focus our analysis on the 
prototypical 1 : 1 dimer chain with beads of different 
masses (denoted mi and m-z). The rescaled equations of 
motion (without gravity; see the discussion below) can 
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be written [l[ 

mi Uj = ( Wj - Uj) k - (Uj - Wj-if , (2) 

m.2Wj — (uj+i — Wj) k — (wj — uj) k , (3) 

where Uj (uij) denotes the displacement of the jth sphere 
of mass mi (m 2 ). We consider a general power-law in- 
teraction to illustrate the comprehensiveness of our ap- 
proach and use Hertzian contact (k = 3/2) to compare 
our theoretical analysis with our numerical and experi- 
mental results. 

The distance from Uj to Wj (from Uj to Uj+i) is D 
(2D). Using D as a small parameter, we develop a long- 
wavelength approximation (LWA) by Taylor-expanding 
Eqs. (O-©, for which we express Uj +1 (lUj-i) as a func- 
tion of uj (wj ) . The resulting PDEs need to be "homog- 
enized" between the two "species" . To accomplish this, 
we follow [lij], postulating a "consistency condition" be- 
tween the two fields: 

w = A (u + b\Du x + b 2 D 2 u xx + b 3 D 3 u xxx 

+b 4 D A u ix + ...). (4) 

We then self-consistently determine the coefficients A and 
bi by demanding that Eqs. ([2|) and ([3]) are identical at 
each order. (Note that the subscripts in the LWA denote 
derivatives.) The parameter A can take the values 1 (for 
acoustic excitations) or -roi/ma (for optical ones). The 
nature of our experimental initial conditions generically 
leads to in-phase waveforms, so we restrict our consid- 
erations to A = 1 hereafter [3^]. Considering the en- 
suing PDE at orders D k -D k+3 , we find that b\ = 1, 
b 2 = rni/(rni +m 2 ), b 3 = (2mi - m 2 )/[3(mi +m 2 )], and 
64 = mi(m\ — m\rti2 + m 2 )/[3(mi + m 2 ) 3 ]. Observe that 
these results correctly capture the uniform-chain limit of 
mi = m 2 , for which w should represent the Taylor ex- 
pansion of radius D around u. 
The resulting PDE, 

U TT - - U x U xx -f - Gu x ^^X ^XX^XXX 

+ iu k r 1 u ix , (5) 

has the same form (with different coefficients) as that 
obtained for uniform chains In Eq. ©, r = 

t^/2kD k+1 / (mi + m 2 ) is a rescaled time, G — D 2 (2 — 
3k + fc 2 )m 2 /[6(mi + m 2 ) 2 ], H = 2D 2 (k - l)(2m\ + 
TOim 2 — mf,)/[6(mi + 7712) ], and / = 2D 2 (m\ — mim 2 + 
m 2 )/[6(mi + ma) 2 ]. We seek traveling-wave solutions 
u = u(£), with £ = x - V s t (and V s = gjp), and ob- 
tain an ordinary differential equation (ODE) for u% — 
v. We then change variables with the transformation 
v = z p , where the power p is chosen so that terms 
proportional to z p ~ 3 z| disappear in the resulting ODE 
for z — z(£). Finally, an integrating factor z a , with 
a = 1 — kp + 3(p — 1) + pH/I, converts the ODE to a 
tractable form, ztt = fiz 11 — az where /i = V 2 /(I(p + a)), 



n = 1 + p — kp, and a — l/[I(kp + a)]. Direct integration 
then yields: 

u ? = v = z p = B cos^ ((3S,) , (6) 

with B = (^/[I3 2 s(s - l)]) 1/(fc_1) , (3 = ^(1 - r?)/2, and 
s = 2/(1-77). 

The existence of such a trigonometric solution in the 
nonlinearly dispersive LWA suggests the possibility of 
finite-width solutions that essentially consist of a single 
arch of the profile of Eq. ©, similar to what has been 
derived for uniform chains [l|. Among the most notable 
properties of such solutions that are testable both experi- 
mentally and numerically are the amplitude- velocity scal- 
ing B ~ Vs^ k (which is similar to the single-species 
case (HQ) an d the solution width n/f3 (which explicitly 
depends on the mass ratio) . We have obtained an analyt- 
ical expression for the width that is valid for all k, but the 
formula is too lengthy to write explicitly in the general 
case. The expression for (3 in the case of Hertzian interac- 
tions (i. e., k = 3/2) is P = (30-8a; + 8a; 2 )- 1 (\/3(8~5tJ + 
5w 2 -V4-4a; + 13w 2 - 18w 3 + 9w 4 )(((l+w) 2 (15-4w + 
4a; 2 )) /(£> 2 (34- 42u; + 59lo 2 - Mto 3 + 17w 4 + (-8 + 5lu - 
5w 2 )V4 - 4w + 13w 2 - 18u 3 + 9a- 4 ))) 1 / 2 ), where u = 
m 2 /mi. This result, which is one of the main findings of 
the present work, generalizes all previously known limit- 
ing cases - namely mi = m 2 , for which (3 — y/lO/(5D) 
[resulting in pulses extending to 57r/VT0 ~ 5 sites], and 
mi ^> m 2 [resulting in pulses of vTOtt ~ 10 sites (com- 
posed of 5 cells with 2 sites each)] 

Comparison between experiments, numerical simula- 
tions, and theory. To begin, we examine the scaling of 
the maximum dynamic force F m of the solitary waves 
versus their propagation speed V s . The theory predicts 
that F m - B k - Vs k/(k ^ 1 \ which for k = 3/2 yields 
V s ~ Fm 6 ■ We tested this experimentally using a 1:1 
steekPTFE dimer chain by dropping the striker (a stain- 
less steel bead) from initial heights ranging from m to 
1.2 m. We placed sensors in beads 15 and 34 and mea- 
sured the corresponding peak force. We averaged their 
amplitudes [F m = (F i5 + F 3 4)/2] and obtained the cor- 
responding wave speed (V s ) using time-of-flight measure- 
ments. We calculated these diagnostics similarly in our 
numerical simulations. As shown in Fig.[3Ja), we obtain 
very good agreement between numerics and experiments. 
A least-squares fit of the numerical simulations (using 
the experimental configuration of 38 beads and removing 
gravity) yields V s ~ F^ 1666 , in excellent agreement with 
the theory (despite the small number of particles). 

To better highlight some of the interesting dynam- 
ics of dimer chains, we examined the evolution of the 
solitary-wave width (the full width at half maximum, 
or FWHM) as the wave progresses down the chain of 
beads. As shown in Fig. [3jb), the numerical and ex- 
perimental results are in very good agreement with each 
other and with the "homogenized" theory (which should 
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FIG. 3: (Color online) Comparison of experiments, numerical 
simulations, and theory, (a) Scaling of the maximum dy- 
namic force F m versus propagation velocity Vi in a 1:1 stain- 
less steekPTFE chain. The numerical simulations are shown 
by the solid curve and the experimental results are shown by 
points (with error bars). The dashed curve shows the numer- 
ical results using the elastic modulus E = 0.6 GPa for PTFE, 
the nominal static value reported in most of the literature 
[iH , [33|]. (b) Evolution of solitary wave width (full width at 
half maximum, or FWHM) as a function of bead number (also 
shown as a log-log plot in the inset). The experimental values 
are shown by solid (green) squares and the numerical values 
are shown by open (red) circles. (In both cases, we include 
curves between the points as visual guides.) The theoretical 
value for the FWHM with mi ^> 7712 is given by line (1), that 
for the 1:1 steebPTFE chain is given by (2), and that for a 
homogeneous chain is given by (3). 



be expected to capture the average of the relevant os- 
cillations). The F m values in the numerics and experi- 
ments alternate from one bead to another because con- 
secutive beads are composed of different materials. As 
one can see from Fig. [3Jb), the average FWHM of the 
dimer is decidely different from that of a homogeneous 
chain (mi = mz) and even from that of the ra\ 3> m,2 
limiting case. It is best captured by the theoretical cal- 
culation above (for the proper ratio of masses), clearly 
evincing the relevance of our theoretical approach. As a 
quantitative measure, the relative error of the theoreti- 
cal prediction versus the computational-average FWHM 
is 22.8% for m x = m 2 , 9.2% for mi > m 2 , and 2.5% for 
the relevant mass ratio. In conjunction with the theo- 
retical analysis developed in this paper, our experiments 
and numerics reveal that the equilibrium conditions (in a 
horizontal chain) for the wave width are achieved even at 
very short propagation distances (i.e., with a very small 
number of beads). For vertical chains of beads that con- 



tain a large number of particles (several hundred or more) 
or are excited by smaller-amplitude pulses, the presence 
of the nonuniform gravitational precomprcssion should 
be taken into account. As discussed in [3], this leads to 
a 1/3 power law scaling in the wave width as a function 
of particle number. 

Conclusions. We examined the propagation of solitary 
waves in heterogeneous, periodic chains of granular me- 
dia using experiments, numerical simulations, and the- 
ory. Using different periodicities, we found that such het- 
erogeneous systems robustly support the formation and 
propagation of highly localized nonlinear solitary waves 
(with appropriate widths that depend on the periodic- 
ity). We used force- velocity scaling (which is the same as 
for homogeneous chains) and solitary-wave width (which 
depends on the mass ratio of the dimer materials) as 
relevant benchmarks for the very good agreement be- 
tween our three approaches. This qualitative and quan- 
titative understanding of the dimer dynamics also paves 
the way for studies in increasingly heterogeneous media 
(e.g., trimers) in both one and higher dimensions, as well 
as of more complex (optical) modes in such systems. 

Acknowledgements We acknowledge support from the 
Caltech Information Science and Technology initiative 
(M.A.P.), start-up funds from Caltech (CD.), and NSF- 
DMS and CAREER (P.G.K.). I.S. performed his work 
as a visiting student of GALCIT at Caltech. We thank 
A. Molinari and V. F. Nesterenko for useful discussions. 

* Electronic address: daraio@caltech.edu 



[1] V. F. Nesterenko, Dynamics of Heterogeneous Materials 

(Springer- Verlag, New York, NY, 2001). 
[2] C. Daraio et al., Phys. Rev. E 73, 026610 (2006). 
[3] C. Coste et al., Phys. Rev. E 56, 6104 (1997). 
[4] E. J. Hinch and S. Saint- Jean, Proc. Royal Soc. London 

A 455, 3201 (1999). 
[5] E. Hascoet and H. J. Herrmann, Eur. Phys. Journal B 

14, 183 (2000). 
[6] M. Manciu et al., Physica D 157, 226 (2001). 
[7] J. Hong and A. Xu, Phys. Rev. E 63, 061310 (2001). 
[8] J. Hong and A. Xu, App. Phys. Lett. 81, 4868 (2002). 
[9] J. Hong, Phys. Rev. Lett. 94, 108001 (2005). 
[10] C. Daraio et al., Phys. Rev. Lett. 96, 058002 (2006). 
[11] R. Doney and S. Sen, Phys. Rev. Lett. 97, 155502 (2006). 
[12] L. Vergara, Phys. Rev. E 73, 066623 (2006). 
[13] C. Daraio et al., Phys. Rev. E 72, 016603 (2005). 
[14] V. F. Nesterenko et al., Phys. Rev. Lett. 95, 158702 
(2005). 

[15] S. Job et al., Phys. Rev. Lett. 94, 178002 (2005). 
[16] A. Sokolow et al., Europhys. Lett. 77 (2007). 
[17] P. Rosenau and J. M. Hyman, Phys. Rev. Lett. 70, 564 
(1993). 

[18] R. Reigada et al., Chaos 13, 646 (2003). 

[19] D. K. Campbell et al., Chaos 15, 015101 (2005). 

[20] J. Ford, Phys. Rep. 213, 271 (1992). 

[21] Y. S. Kivshar and G. P. Agrawal, Optical Solitons: From 



5 



Fibers to Photonic Crystals (Academic Press, San Diego, 

California, 2003). 
[22] J. Fleischer et al., Opt. Express 13, 1780 (2005). 
[23] O. Morsch and M. Oberthaler, Rev. Mod. Phys. 78, 179 

(2006). 

[24] M. Peyrard, Nonlinearity 17, Rl (2004). 

[25] H. Bilz et al., Phys. Rev. Lett. 48, 264 (1982). 

[26] P. C. Dash and K. Patnaik, Prog. Theor. Phys. 65, 1526 

(1981) . 

[27] M. J. Rice and E. J. Mele, Phys. Rev. Lett. 49, 1455 

(1982) . 

[28] A. A. Sukhorukov and Y. S. Kivshar, Opt. Lett. 27, 2112 
(2002); R. Morandotti et al., Opt. Lett. 29, 2890 (2002). 
[29] M. Sato et al., Rev. Mod. Phys. 78, 137 (2006). 



[30] A.-C. Hladky-Hennion and G. Allan, J. of App. Phys. 
98, 054909 (2005). 

[31] ASM Metals Reference Book (American Society for Met- 
als, Metals Park, OH, 1983), 2nd ed. 

[32] http://www.efunda.com] 

[33] www . dupont . com/teflon/chemical/ 

[34] S. Pnevmatikos et al., Phys. Rev. B 33, 2308 (1986). 

[35] To investigate optical excitations, it would be necessary 
to employ a fundamentally different and more complex 
framework involving a single precompressed bead and 
excitations of prescribed frequency in the corresponding 
linear-spectrum gap of such a setting. 



